Antiferromagnetic spin excitations in single crystals of nonsuperconducting Lii^FeAs 
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We use neutron scattering to determine spin excitations in single crystals of nonsuperconducting 
Lii-zFeAs throughout the Brillouin zone. Although angle resolved photoemission experiments and 
local density approximation calculations suggest poor Fermi surface nesting conditions for antifer- 
romagnetic (AF) order, spin excitations in Lii-^FeAs occur at the AF wave vectors Q = (1,0) 
at low energies, but move to wave vectors Q = (±0.5, ±0.5) near the zone boundary with a total 
magnetic bandwidth comparable to that of BaFe2As2- These results reveal that AF spin excitations 
still dominate the low-energy physics of these materials and suggest both itinerancy and strong 
electron-electron correlations are essential to understand the measured magnetic excitations. 



PACS numbers: 74.25.Ha, 74.70.-b, 78.70.Nx 

Understanding whether magnetism is responsible for 
superconductivity in FeAs-based materials continues to 
be one of the most important unresolved problems in 
modern condensed matter physics p~H3]. For a typical 
iron arsenide such as LaFeAsO [1] , band structure calcu- 
lations predict the presence of the hole-like Fermi surfaces 
at the r(0, 0) point and electron-like Fermi surfaces at the 
M(l, 0)/(0, 1) points in the Brillouin zone [Fig. 1(a)] [SJ. 
As a consequence, Fermi surface nesting and quasiparti- 
cle excitations between the hole and electron pockets can 
give rise to static antiferromagnetic (AF) spin-density- 
wave order at the in-plane wave vector Q — (1,0) [TJ. 
Indeed, neutron scattering experiments have shown the 
presence of the Q = (1,0) AF order in the parent com- 
pounds of iron arsenide superconductors, and doping to 
induce superconductivity suppresses the static AF order 
[5]. In addition, angle resolved photoemission measure- 
ments [5] have confirmed the expected hole and electron 
pockets in superconducting iron arsenides, thus providing 
evidence for superconductivity arising from the sign re- 
vised electron-hole inter-pocket quasiparticle excitations 

EH EM. 

Of all the FeAs-based superconductors pQ, LiFeAs is 
special since it has the highest transition temperature 
(T c = 18 K) amongst the stoichiometric compounds [T3T - 
H7| . Furthermore, it does not have static AF order due 
to the poor Fermi surface nesting properties with shallow 
hole pockets near the T(0,0) [IS]. It has been suggested 
that the flat tops of the hole pockets in LiFeAs imply 
a large density of states near the Fermi surface, which 
should promote ferromagnetic (FM) , instead of the usual 
AF, spin fluctuations for superconductivity |19| . If this 
is indeed the case, AF spin fluctuations should not be 
fundamental to the superconductivity of FeAs-based ma- 



terials and the superconducting pairing would not be in 
the spin singlet channel. A determination of the mag- 
netic properties in LiFeAs is thus important to complete 
our understanding about the role of magnetism in the 
superconductivity of FeAs-based materials. 

In this paper, we present inelastic neutron scattering 
measurements on single crystals of nonsuperconducting 
Lii_ x FeAs with x = 0.06 ± 0.01, where there is no static 
AF order. As a function of increasing energy, spin ex- 
citations in Lio^FeAs have a spin gap below A ss 13 
meV, are centered at the AF wave vector Q = (1,0) 
for energies up to ~80 meV, and then split into two 
vertical bands of scattering before moving into the zone 
boundaries at the wave vectors Q' = (±1/2, ±1/2) near 
E ss 130 meV. These Q' vectors have been observed in 
the spin excitations of FeTe/Se compounds and imply 
the existence of strong competition between FM and AF 
exchange couplings [20]. While the dispersions of the 
low-energy spin excitations (E < 80 meV) in Li .94FeAs 
are similar to that of (Ba,Ca,Sr)Fe 2 As 2 Q] [211 [23], the 
high-energy spin excitations near the zone boundary are 
quite different from these materials, and cannot be mod- 
eled from a simple Heisenberg Hamiltonian with effective 
nearest (J± a and Jib) and next nearest neighbor (J 2 ) ex- 
change couplings [2 [3T]. By integrating the local sus- 
ceptibility x"(w) in absolute units over the entire band- 
width of spin excitations, we find the spin fluctuating 
moment (m 2 ) = 2.1 ± 0.6 n 2 B , a value that is compara- 
ble with other pnictides. Therefore, spin excitations in 
Li .94FeAs are similar to other iron pnictides but are not 
directly associated with Fermi surface nesting from hole 
and electron pockets, contrary to expectations from local 
density approximation calculations [TSJ US] • 

Our experiments were carried out on the ARCS time- 
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FIG. 1: (color online), (a) Fermi surfaces from the spin- 
restricted local density approximation calculation for LiFeAs 
|18l I19j . There are two hole-like Fermi surfaces near the 
F(0,0) point with d yz /d X z character and one electron-like 
Fermi surface near the M(1,0) point. The nesting condi- 
tion for the expected AF nesting wave vector Qafm = (1,0) 
is not favorable [18] . (b) Zero field cooled (ZFC) and field 
cooled (FC) susceptibility measurements on Lio.94FeAs. No 
superconductivity was observed due to Li deficiency, (c) The 
dashed lines show spin wave dispersions along the [H, 0] and 
[1,K] directions for BaFe 2 As 2 at 5 K [21 . The filled cir- 
cles show the measured spin excitation dispersions along the 
[H, 0] and [1,K] directions for Lio.94FeAs. While spin waves 
in BaFe2As2 extend up to 200 meV along the [1, K] direction, 
spin excitations in Lio.g4FeAs reach the zone boundary near 
Q= (1,0.5). 



of-flight chopper spectrometer at the spallation neutron 
source, Oak Ridge National laboratory. We also per- 
formed thermal triple-axis spectrometer measurements 
on the BT-7 triple-axis spectrometer at NIST center for 
neutron Research. Our single crystals were grown using 
the flux method and inductively coupled plasma analy- 
sis on the samples showed that the compositions of the 
crystals are Li .94±o.oiFeAs. Figure 1(b) shows zero field 
cooled (ZFC) and field cooled (FC) susceptibility mea- 
surements on Lio.g4FcAs, which indicate spin glass be- 
havior with no evidence for superconductivity. To study 
the spin excitations, we co-aligned 7.5 g of single crystals 
of Lio.94FeAs (with a mosaic of 2°) and loaded the sam- 
ples inside a He refrigerator or cryostat. To facilitate easy 
comparison with spin wave measurements in BaFe2As2 
[21j . we denned the wave vector Q at (q x , q y , q z ) as 
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FIG. 2: (color online). Triple-axis measurements to search for 
static AF order and spin excitations in Lio.94FeAs. (a) Elastic 
scattering along the [H, 0, 3] and (b) [1, 0, L] directions at 2 K 
show no evidence of AF order at the expected position Q = 
(1, 0, 3). The arrows indicate Al sample holder scattering, (c) 
Constant-Q scans at the wave vectors Q = (1,0,3) (signal) 
and Q = (0.4, 0, 3) (background) positions at 2 K. A clear 
spin gap is seen at A = 13 meV. (d) Constant-energy scans at 
E = 9, 16 meV along the [H, 0, 3] direction. While the scan at 
E — 9 meV is featureless, a clear peak is seen at E = 16 meV 
confirming the spin gap. (e) Imaginary part of the dynamic 
susceptibility x" at 2 K and 190 K. The magnitude of the 
spin gap is unchanged between 2 and 190 K. (f) Temperature 
dependence of the x"(Q) at E = 16 meV for 2, 100, and 190 
K. The x"(Q) is almost temperature independent between 2 K 
and 190 K. Error bars where indicated represent one standard 
deviation. 



(H,K,L) = (q x a/2Tr,q y b/2ir,q z c/2Tr) reciprocal lattice 
units (rlu), where a = b = 5.316 A, and c = 6.306 A. 
For both triple-axis and ARCS measurements, we aligned 
crystals in the [H, 0, L] scattering zone. The ARCS data 
are normalized to absolute units using a vanadium stan- 
dard. The incident beam energies were = 80, 140, 250 
meV with Ei parallel to the c-axis. 

Before describing in detail the spin excitation disper- 
sion curves and dynamic local susceptibility in Figs. 1 
(c)-l(e), we first discuss the triple-axis measurements on 
the static AF order and spin excitations. Figures 1(a) 
and 1(b) show elastic scattering along the [H, 0, 3] and 
[1,0, L] directions across the expected AF peak position 
(1,0,3), respectively. In contrast to Nai^FeAs, where 
static AF order is clearly observed [2"3] . there is no evi- 
dence for static AF order in this sample. To search for 
AF spin excitations, we carried out constant-Q scans 
at the AF wave vector Q — (1,0,3) and background 
(0.4,0,3) positions. The outcome in Fig. 2(c) shows 
a step-like increase in scattering above background for 
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FIG. 3: (color online). Constant-energy images of the scat- 
tering in the [H, K] zone as a function of increasing energy 
for Lio.94FeAs at energy transfers of (a) E = 25 ± 5 meV; 
(b) 45 ± 5 meV (with E, = 80 meV); (c) 70 ± 10 mV; (d) 
90±10meV; (e) llOilOmeV; (f) 130±10meV; (g) 150 ±10 
meV; (h) 170 ± 10 meV, all with 2% = 250 meV. The scatter- 
ing intensity is in absolute units. The box in (b) shows the 
Brillouin zone used to integrated the susceptibility. 



E > 13 meV, clearly suggesting the presence of a large 
spin gap of A = 13 meV. To confirm there is indeed 
a spin gap, we carried out constant-energy scans along 
the [H, 0, 3] direction at E = 9 and 16 meV as shown in 
Fig. 1(d). While the scattering is featureless at E = 9 
meV, there is a clear peak centered at Q = (1,0,3) at 
E = 16 meV. Figure 3(e) shows temperature dependence 
of the imaginary part of dynamic susceptibility x" {E) ob- 
tained by subtracting the background and correcting for 
the Bose population factor. Surprisingly, the spin gap has 
no observable temperature dependence between 2 K and 
190 K, much different from the temperature dependence 
of the spin gaps in the (Ba,Sr,Ca)Fe2As2 family of ma- 
terials |25H27j . which disappear rapidly with increasing 
temperature. The weak temperature dependence of the 
dynamic susceptibility has been confirmed by constant- 
energy scans in Fig. 2(f), where x"(Q) at -E 1 = 16 meV 
remains essentially unchanged from 2 K to 190 K. 

Figure 3 summarizes the ARCS time-of-flight mea- 
surements on Lio.94FeAs at 5 K. Since spin excitations 
in Lio.94FeAs have no c-axis modulations, we show in 
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FIG. 4: (color online). Constant-energy cuts of the spin exci- 
tation dispersion as a function of increasing energy along the 
[H,0] and [1,K] directions for Lio.94FeAs. The dashed curves 
show identical cuts for spin waves of BaFe2As2 divided by 2 
[21] , Both are in absolute units. Constant-energy cuts along 
the [H, 0] direction at (a) 55 ± 5, (b) 75 ± 5, (c) 95 ± 10, 
and (d) 135 ± 10 meV. Similar cuts along the [1, K, ] direction 
are shown in (e)-(h). The dynamic spin correlation length is 
£^ 12 ±3 A. 



Figs. 3(a)-3(h) two-dimensional constant-energy (E) im- 
ages of the scattering in the (H, K) plane for E = 25 ± 5, 
45 ± 5, 70 ± 10, 90 ± 10, 110 ± 10, 130 ± 10, 150 ± 10, 
and 170 ± 10 meV, respectively. For energies between 
25 ± 5 < E < 90 ± 10 meV, spin excitations form trans- 
versely elongated ellipses centered around AF Q = (1,0). 
The intensity of spin excitations decrease with increas- 
ing energy, this is remarkably similar to spin waves in 
BaFe2As2 [2T]. For energies above 90 meV, spin excita- 
tions split into two horizontal arcs that separate further 
with increasing energy. The excitations finally merge into 
Q = (m±0.5, n±0.5) (m, n = 0, 1, 2) and become weaker 
above 150 meV in Fig. 3(g) and 3(h). 

In order to determine the dispersion of spin excita- 
tions for Lio.94 FeAs, we show in Fig. 4 cuts through 
the two-dimensional images in Fig. 3 and compare with 
identical cuts for BaFe2As2- Figures 4(a)-4(d) show 
constant-energy cuts along the [H, 0] direction for ener- 
gies of E = 55 ± 5, 75 ± 5, 95 ± 5, 135 ± 5 meV, respec- 
tively, while the dashed lines show identical spin wave 
cuts for BaFe2As2 [H]. Since both measurements were 
taken in absolute units, we can see that spin excitations 
in Lio.94 FeAs are similar to that of BaFe2As2 per Fe be- 
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low 95 meV [3T]. Figures 4(e)-4(h) show constant-energy 
cuts along the [1,K] direction for identical energies as 
that of Figs. 4(a)-4(d). For energies above 95 meV, the 
strength of the spin excitations in Lig^FeAs are rapidly 
suppressed compared to those of BaFe2As2 and become 
very weak above E — 135 meV. This can originate from 
an absence of magnetic scattering, or that the scattering 
is very broad as might occur when an itinerant electron 
system interacting with Stoner excitations. This is dif- 
ferent from spin waves in BaFe2As2, which extend up to 
250 meV. Based on these constant-energy cuts, we show 
in Figs. 1(c) and 1(d) the comparison of spin excitation 
dispersions of Lio.g4FeAs (filled circles) with those of spin 
waves in BaFe2As2 (dashed lines). They are similar for 
energies between 50 — 95 meV, while the spin excitations 
in Lio.94FeAs are broader below 50 meV. 

We have attempted, but failed, to use a simple Heisen- 
berg Hamiltonian with the effective nearest and next 
nearest neighbor exchange couplings to fit the observed 
spin excitations spectra [1] [21] . For all the possible com- 
binations of the J i a , Ju,, and J 2 , the expected zone 
boundary spin excitations are quite different from the 
observed spectra (see supplementary information) [2"5] , 
If we include the next-next nearest neighbor exchange 
coupling J3, the expected spectra near the zone bound- 
ary have some resemblance to the data in Fig. 3 al- 
though the low-energy excitations would be different (see 
supplementary information). This means that the ef- 
fective exchange couplings in Lio.94FeAs are extremely 
long-ranged, a hallmark that itinerant electrons are im- 
portant for spin excitations in this material. Since the 
data close to the band top along the [1,K] direction are 
higher in energy than along the [H, 0] direction, we need 
Jib < to recover this feature in a Ji a -J\b-J2-J-i model 
(see supplementary information). This means that effec- 
tive exchange interactions in Lio.94FeAs may be similar 
to the (Ca,Sr,Ba)Fe 2 As 2 iron pnictides [TJ[STJ[33] in spite 
of their different zone boundary spectra. 

Finally, we show in Figure 1(e) the energy depen- 
dence of the local susceptibility, defined as x' '(E) — 
J ~x!'((\,E)d(\/ J dq, where the average is over the mag- 
netic scattering signal x" (l> E) over the Brillouin zone 
[Fig. 3(b)] |29j. The corresponding fluctuating moment 
(m 2 ^ = 2.1 ± 0.6 ii 2 B per formula unit. We can use 
both pure local and itinerant spin models to sketch a 
basic physical picture based on the moment value. If 
we assume a quantum local spin model to describe the 
fluctuations, the moment value implies the spin value is 
about one. If we take a pure itinerant model, our result 
would suggest that at least three electrons per iron site 
occupy the states with energies up to the magnetic band- 
width (~ 150 meV) below the Fermi energy. This sug- 
gests that the bandwidths of the electron bands near the 
Fermi surface are extremely narrow. In other words, the 
band renormalization factors are large and the electron- 
electron correlations must be strong. 



In summary, we measured spin excitations in single 
crystals of Lio.g4FeAs. Similar to other iron pnictides, the 
low energy excitations are still strongly AF [30]. How- 
ever, comparing to other iron pnictides, they have sev- 
eral distinct properties: (a) a larger spin gap, close to 
13 meV that is essentially temperature independent be- 
low 190 K; (b) a comparable total magnetic bandwidth; 
(c) different wave vectors at the zone boundary for high 
energy excitations. Moreover, the excitations can not 
be described by magnetic models with only short range 
magnetic exchange couplings. Our results suggest the 
AF spin fluctuations are fundamental to the supercon- 
ductivity of FeAs-based materials. FM fluctuations exist 
in Lio.94FeAs, but they only affect the high energy spin 
excitations. 

During the process of writing up this paper, we became 
aware of a related work on powder samples of supercon- 
ducting LiFeAs, where AF spin fluctuations have been 
reported [3"T] . 
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netic spin excitations in single crystals of 
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In order to test if a simple Jxa-Jib-J^-Js Heisenberg 
Hamiltonian can reproduce the observed spin excitations 
spectra in Fig. 3, we simulated the expected spin wave 
spectra using Heisenberg Hamiltonian [1J. To facilitate 
direct comparison with the data in Fig. 3, we normal- 
ized the calculated intesity at 90 meV to be the same as 
that in Fig. 3. Therefore, the calculated spectra can be 
directly compared with the observed spectra. Figure SI5 
shows spin wave calculations assuming SJ\ a = 23 meV, 
SJ lb = -4 meV, SJ 2 = 13 meV, and SJ 3 = -5 meV. 
While the zone boundary spectra have some similarity 
to the data, the spectra clearly disagree with the data 
around intermediate energies. Figure SI6 shows similar 
calculation assuming SJ\ a = SJ\b = 10 meV, SJ 2 = 20 
meV, and SJ3 = —8 meV. Figure SI7 plots calculations 
assuming SJi a = 23 meV, S Jib = —4 meV, SJ 2 = 13 



In the model, the dispersion along [H, 0] -direction sees 
the maximum at (1/2, 0), while along [1, iT]-direction, the 
maximum is at (l,q) where q ~ 1/2 is in fact parameter 
dependent. But one should not only compare w(l/2, 0) 
and oj(q, 0) to find which direction reaches a higher top, 
because there is twinning. Once we have twinning into 
play, we also need to compare the band near (0, 1). Mark 
that in discussing this region, H and K directions should 
be interchanged when compared with the experimental 
frame. Starting from (0, 1), the dispersion reaches maxi- 
mum along if -direction at (0, 1/2) and along H direction 
at (1/2, 1). Therefore, in order to see what we see in the 



meV, and SJ3 = meV; and Figure SI8 show calcula- 
tion with SJ\ a = SJu = 10 meV, SJ 2 = 20 meV, and 
SJ3 — meV. The spin wave band tops are in Figs. 
SI5-SI8 are 150 meV, 150 meV, 110 meV, and 100 meV, 
respectively. To obtain similar scattering pattern as ob- 
served near the band top by a Heisenberg Hamiltonian 
with effective exchange couplings, the next-next nearest 
neighbor exchange coupling J 3 has to be included. 

We note that spin excitations close to the band top 
in Figl(e) clearly shows that zone boundary energy in 
the [1, K] direction is higher in energy than that along 
the [H, 0] direction. In the following we show that in a 
■j ha~ hb- h~ h model with antiferromagnetic Ji a and J 2 
and ferromagnetic J 3 , we need and only need Jib < 
(ferromagnetic) to recover this feature. 

The Ji a -J\b-J2-J-i model is analogous to the Ji a -J\b~J2 
spin model used to describe CaFe2 AS2 [1] . In the following 
we will consider a detwinned system, and one should bear 
in mind that one q-point (q x ,q y ) in a twinned system 
corresponds to (q Xl q y ) and (q y ,q x ) in a detwinned one. 
The spinwave dispersion for this model is given by 

u(h, k) = Sy,4 2 (M)- B 2 (h,k), (1) 

where 



(2) 
I 

experiment, we should have 

max[w(l/2,0),u;(0,l/2)] < max[w(l, q), w(l/2, 1)]. (3) 

Now we make a statement and prove it: < is 
sufficient and necessary for the Eq.(|3]) to hold. First we 
prove the sufficiency. 

w 2 (l/2, 1) - w 2 (l/2, 0) = 16J 16 (- J la + J u - 2J 2 + 2J 3 ).(4) 

From this we know if J^ < 0, we have u>(l/2, 1) > 
w(l/2,0). On the other hand, 



A(h,k) = 2[Ji 6 (cos(ttA;) - 1) + Ji„ + 2Ji + J 3 (cos(27r/i) + cos(27rfc) - 2)], 
B(h,k) = 2( Ji a cos(7r/i) + 2J 2 cos(nh) cos(7rfc)). 



w 2 (l/2, 1) - uj 2 (0, 1/2) - 4(J 2 Q - 2J la J u + Jib(3J u - 4J 2 + 4J 3 )). (5) 

From this we know if 

Jib < + 2J 2 - 2J 3 - ^{J la + 2J 2 -2J 3 y-ZJl a , (6) 

cj(1/2,1) > w(Q,l/2). 
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But of course 

~(J U + 2J 2 - 2J 3 - y(Ji + 2J 2 -2J 3 ) 2 -3J 1 2 a > 0,(7) 
therefore Jib < is sufficient to make the highest energy 



along [l,K] direction higher than [H,0] direction. 

Then we prove the necessity, ft is a proof by contra- 
diction. Suppose Jib > 0, then from above we know that 
ui(l/2, 1) < w(l/2,0). Also notice that when J lb > 0, 



oj 2 (1 , q y ) = 4(J lo + 2J 2 + J 16 (cos(7rg) - f ) - 2J 3 sm 2 (7rg)) 2 - 4(J lo + 2J 2 cos(7rq)) 2 (8) 

< 4(J lo + 2J 2 - 2J 3 ) 2 - 4( J la - 2 J 2 ) 2 
= W(Ji a - J 3 )(2J 2 - J 3 ) 

< 4( J lQ + 2J 2 - 2J 3 ) 2 
= ^ 2 (V2,0). 




Therefore 

max[w(l/2,0),w(0,l/2)] > max[w(l, g), w(l/2, 1)]. (9) 

Therefore we have proved if Jib > then the highest 
energy along [1,-K] direction is lower than the highest 
energy along [H, 0] direction, i.e., Jib < is necessary to 
recover the feature observed in experimental data. 



f = 45 meV E = 70ineV 




f = 45 meV 




-2 2 -2 2 



[H, 0] (r.l.u) [1,K](r.l.u) 

SI 6: (color online). Constant-energy images of the scatter- 
ing in the [H, K] zone as a function of increasing energy for 
Lio.94FeAs at energy transfers of (a) 45 meV; (b) 70 meV; (c) 
90 mV; (d) 110 meV; (e) 130 meV; (f) 150 meV; (g,h) The 
expected spin wave dispersion along the [H, 0] and [l,K] di- 
rections for LiFeAs. The exchange couplings used to obtain 
these imagines are SJi a = 10 meV, SJib ~ 10 meV, SJ2 = 20 
meV, SJ3 — —8 meV, and anisotropy factor SJ S = 0.08 meV. 





SI 5: (color online). Constant-energy images of the scatter- 
ing in the [H, K] zone as a function of increasing energy for 
Lio.94\FeAs at energy transfers of (a) 45 meV; (b) 70 meV; (c) 
90 mV; (d) 110 meV; (e) 130 meV; (f) 150 meV; (g,h) The 
expected spin wave dispersion along the [H, 0] and [1,K] di- 
rections for LiFeAs. The exchange couplings used to obtain 

these imagines are S J\ a = 23 meV, S Jif, = —4 meV, SJ2 = 13 
meV, SJ3 — —5 meV, and anisotropy factor SJ S = 0.08 meV. 

Spin wave damping is assumed to be F = 0.15E throughout * Electronic address: pdai@ utk.edu| 

the supplementary information. I 1 ] J - zh ao et ai, Nat. Phys. 5, 555 (2009). 
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SI 7: (color online). Constant-energy images of the scatter- 
ing in the [H, K] zone as a function of increasing energy for 
Lio.94FeAs at energy transfers of (a) 45 meV; (b) 70 meV; (c) 
90 mV; (d) 110 meV; (e,f) The expected spin wave disper- 
sion along the [H, 0] and [1, K] directions for LiFeAs. The ex- 
change couplings used to obtain these imagines are S Ji a = 23 
meV, SJ lb = -4 meV, SJ 2 = 13 meV, SJ 3 = meV, and 
anisotropy factor SJ 3 = 0.08 meV. 




-2 

[H, 0] (r.l.u) 



-2 

[1,K] (r.l.u) 



SI 8: (color online). Constant-energy images of the scatter- 
ing in the [H, K] zone as a function of increasing energy for 
Lio.94FeAs at energy transfers of (a) 45 meV; (b) 70 meV; (c) 
90 mV; (d) 100 meV; (e,f) The expected spin wave disper- 
sion along the [H, 0] and [1, K] directions for LiFeAs. The ex- 
change couplings used to obtain these imagines are S Ji a = 10 
meV, SJ lb = 10 meV, SJ 2 = 20 meV, SJ 3 = meV, and 
anisotropy factor SJ 3 = 0.08 meV. 



